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We present a measurement of the branching ratios hom the 6P3/2 state of Ball into all dipole- 
allowed decay channels (6S1/2, 5D3/2 and 505/2)- Measurements were performed on sing le "«Ba+ 
ions in a linear Paul trap with a frequency-doubled mode-locked Ti:Sapphire laser resonant with the 
6S1/2 ^ 6P3/2 transition at 455 nm by detection of electron shelving into the dark 5D5/2 state. By 
driving a n Rabi rotation with a single femtosecond pulse, an absolute measurement of the branching 
ratio to 5D5/2 state was performed. Combined with a measurement of the relative decay rates into 
5D3/2 and 5D5/2 states performed with long trains of highly attenuated 455 nm pulses, it allowed the 
extraction of the absolute ratios of the other two decays. Relative strengths normalized to unity are 
found to be 0.756±0.046, 0.0290±0.0015 and 0.215±0.0064 for 6S1/2, 5D3/2 and 5D5/2 respectively 
This approximately constitutes a threefold improvement over the best previous measurements and 
is a sufficient level of precision to compare to calculated values for dipole matrix elements. 



PACS numbers: 31.15.bw, 32.70.Cs, 32.70.Fw 

Single trapped ions are a valuable physical system 
for many applications including quantum computation 
[U E], frequency standards, optical metrology [3], pre- 
cision searches for drifts in fundamental constants |1] 
and tests of exotic physical theories [S]. Among the ad- 
vantages over alternatives are their long trapping life- 
times and the relative ease of confining single ions to 
a small volume in a trap, thereby reducing the system- 
atic effects and negating the need for quantum statistics. 
The barium ion, particularly the odd 137 isotope with 
nuclear spin 3/2, has been proposed for use in quan- 
tum computation schemes with the hyperfine levels of 
the 6S1/2 ground state for the qubit, as an optical fre- 
quency standard with a 2051 nm clock transition from 
6S1/2, F = 2, mp = to 5D3/2, F' = 0, and as a test 
of parity-nonconservation with a small dipole coupling 
between the otherwise dipole- for bidden 6S1/2 ^ 5D3/2 
transition [3 El [7j. 

Accurate models of atomic wave functions which in- 
clude many-body interactions are necessary to calculate 
dipole and quadrupole matrix elements that appear in 
the calculations of transition rates, energy level shifts and 
line widths in the experiments mentioned above. Mea- 
surements of branching ratios represent a better quantity 
from which to verify such values than, for example, pre- 
cise measurements of the lifetimes of metastable states. 
They are less prone to systematic uncertainties such as 
background gas quenching and stray fields to which the 
long waiting times (tens of seconds) required to accu- 
rately measure lifetimes are sensitive. Here we present 
a single-ion measurement of the branching ratios from 
the 6P3/2 state of ^•^^Ba"'' to the three states allowed via 
dipole transitions, 6S1/2, 5D3/2 and 5D5/2. 

A schematic of the optical and electronic arrangement 
of the experimental apparatus can be found in Fig. 1. 
The ion trap itself is a linear Paul trap with radio- 



frequency quadrupole confining potential and DC voltage 
end caps in ultra high vacuum with operating pressures 
of about 10^^^ torr. The trap dimensions are ~0.5 mm 
radially and ~ 3.3 mm axially. At 0.5 W of inductively- 
coupled RF power at ~ 32 MHz and 100 V end cap po- 
tential, the trap secular frequencies are measured to be 
approximately cjfl/27r = 5 MHz and uJz/'i'K = 2 MHz. 
A Xenon flash lamp is focussed into the trap to pho- 
toionize a thermal beam of atomic barium. The bar- 
ium ion has a lambda-type energy spectrum (shown in 
Fig. 2) with a strong 6S1/2 ^ 6P1/2 transition at 493.4 
nm which is used for both Doppler cooling and fluo- 
rescence detection, and a repump from the metastable 
5D3/2 state at 649.7 nm. Both cooling and repumping 
are accomplished with commercial external-cavity diode 
lasers (Toptica TA-SHG 110 and DL 100 respectively). 
The 493 nm light is generated by frequency doubling 
the output of a 986 nm diode with a KNb03 crystal in 
a bow-tie enhancement cavity. Frequency stabilization 
is accomplished via Doppler-free spectroscopy on vapor 
cells, Te2 for the 493 nm laser and I2 for the 650 nm 
laser. Both lasers are matched onto a dichroic mirror and 
flber-coupled for mode cleaning and to ensure coUinear- 
ity at the trap. A magnetic field of '--^ 4 G applied par- 
allel to the cooling beams breaks the degeneracy of Zee- 
man sublevels. Effective cooling is achieved with approx- 
imately 100 /iW of blue light and 40 of red. State 
discrimination is accomplished by observing ion fluores- 
cence on the cooling transition collected with a Mitutoyo 
N. A. =0.28 objective with a Hamamatsu PMT behind a 
493 nm interference filter. Using a 100 ms detection in- 
terval, "bright" and "dark" histograms of photon counts 
are typically 30 and 200, respectively. This is a suffi- 
cient separation to completely rule out the possibility of 
incorrect state determination. 

Excitation from the ground state to the 6P3/2 state is 
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accomplished with broadband pulses with central wave- 
length of approximately 455 nm generated by frequency 
doubling a mode-locked Ti:Sapphire laser operating at 
910 nm via single pass of a 5 mm /3BaB204 crystal. 
Pulses had a duration of 400 fs before doubling as mea- 
sured by intensity autocorrelation. Pulse energy was ad- 
justed by a DAC-controUed voltage variable RF attenua- 
tor (Minicircuits ZX73-2500) in line between the switch- 
ing AOM and its driver. For the experiment measur- 
ing the relative branching of the two D levels, a time- 
averaged power level of between 1.6 and 3.4 mW of 
455 nm light was sufficient. Single pulse selection for 
the Rabi rotation experiment was accomplished with the 
same AOM synchronized to the 76 MHz repetition rate of 
the laser using an intracavity-monitoring fast diode cou- 
pled to a Potato Semiconductor PO74G08A, 1.125 GHz 
bandwidth AND-gate, and an HP 8013B TTL Pulse Gen- 
erator. Extinction of neighboring pulses was better than 
97% at 455 nm, and asynchronicity was monitored with 
an HP 5335A universal counter and found to occur an av- 
erage of 3-5% of the time. The maximum achievable sin- 
gle pulse energy was 1.68 nJ after transmission through 
the pulse selector. Shot-to-shot variability of the pulse 
energy was approximately 10% in the blue. De-shelving 
from 5D5/2 was accomplished with 50 ms flashes from a 
614 nm high-power LED. 
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FIG. 1: Schematic overview of the apparatus used in this 
experiment. After stabilization to their respective gas cells, 
cooling (493 nm) and repump lasers (650 nm) are coupled on 
a dichroic mirror and fiber coupled for mode-cleaning and to 
ensure coUinearity. Detection of ion fluorescence, after col- 
lection optics and flltering through a 493 nm narrow band 
interference fllter, can be switched between a cooled electron- 
multiplied CCD (EMCCD) Andor iXon camera with sin- 
gle photon-sensitivity or a Hamamatsu photomultiplier tube 
(PMT). The Xenon flash lamp has high spectral content in 
the UV and thereby photoionizes of a thermal beam of atomic 
barium. The Ti:Sapphire, independently synchronized to the 
measurement apparatus through the 1.125 GHz gate and HP 
8031B pulse generator, provides 400 fs pulses at 910 nm cen- 
tral wavelength which are doubled to 455 nm after a single 
pass through the BBO crystal. The 614 nm high power LED 
de-shelves the ion from the 5D5/2 state. 
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FIG. 2: Relevant energy levels and lifetimes in the ^^^Ba"*" ion. 
D state lifetimes are experimental values [23]. The decays 
under consideration are from 6P3/2 to 6S1/2 at 455.403 nm, 
5D3/2 at 585.530 nm (not shown in digram) and 5D5/2 at 
614.171 nm. 



To establish a value for the relative strengths of the 
decays to the two D states, the ion was prepared in the 
6S1/2 state by optically pumping out of the 5D3/2 state 
with the 650 nm laser. With the two cooling lasers shut- 
tered off, trains of highly attenuated 455 nm pulses of 
various duration from 0.1 to 10 ms and average powers 
of 1.6 to 3.4 mW were focused onto the ion, weakly ex- 
citing it to the 6P3/2 state, which has a characteristic 
lifetime of 6.32 ns [5]. After this state decayed, the cool- 
ing lasers were turned on to look for ion fluorescence. 
The absence of photons at 493 nm indicated that the ion 
was shelved to the 5D5/2 state, but decays to 6S1/2 and 
5D3/2 were indistinguishable using this method. Plotting 
the probability of finding the ion in the dark state as a 
function of the exposure time to 455 nm light, one finds 
that the dark probability asymptotically approaches a 
constant value corresponding to the ratio of strengths of 
decay to the two D states as displayed in Fig. 3. The 
solid line is a fit to the data using the expected form of 
this function, Pdark = A + B (l + e^^*) where A is any 
offset from zero due to false dark counts, mostly resulting 
from collisions of the ion with background gas particles 
and typically much less than 1%, B is the desired ratio 
and A is a parameter which depends on laser intensity 
and proximity to resonance. Using this excitation and 
detection scheme, we derive the relative strength of D 
state decays from this asymptotic value of the dark state 
probabihty as a function of shelving pulse length. Values 
for the ratios excluding decays to the ground state are 
0.881 and 0.119 for 5D5/2 and 5D3/2 respectively. Each 
has a 3% statistical and 0.5% systematic uncertainty, the 
latter arising mainly from leakage of cooling light during 
the shelving phase of the experiment and coUisional de- 
excitation from the 5D5/2 state due to background gas. 
Optical attenuation of the cooling light was at least 85 
dB, which ensures that the systematic decrease in the 
measured ratio is less than 0.5%. From extrapolated 
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quenching rates at our operating pressure of 2.0 x 10""'^"'^ 
torr over the 100 ms duration of each experiment, coUi- 
sional de-excitation contributes to a systematically low 
value by of order 0.1% 0. 



1.0-1 




0.0-1 . , . , . , . , . ,— 

0.000 0.002 0.004 0.006 0.008 0,010 

Exposure time (s) 



FIG. 3: Representative data obtained to establish a value for 
the relative strengths of the decays to the two D states in 
^^*Ba^. The ion was exposed to highly attenuated 455 nm 
light, weakly driving it to the 6P3/2 state. At long exposure 
times, the probability that the ion is dark will determine the 
relative proportion of the trials where the decay from 6P3/2 
was to 5D5/2- The solid line is a fit to the data using a three- 
parameter model, Pdark ~ A -f B (1 -I- e^"^*), where the fit 
parameter B determines the ratio of the relative strengths of 
the D state decays. 

In order to establish absolute values for the decay ra- 
tios to all three states, one must devise a means to inde- 
pendently distinguish the 6S1/2 and 5D3/2 decays. This 
is accomplished by driving the ground state to 6P3/2 
transition with single 400 fs pulses. Similar excitation 
schemes have been employed on the 5Si/2 5Pi/2 tran- 
sition in atomic *^Rb [TD] and on the 5Si/2 5P3/2 
in i"Cd+ and 6S1/2 ^ 6P1/2 in i^iYb+ ions for use 
in coherent coupling of ion spin and photon polariza- 
tion [TT| [12 • The pulses rotate the state of the ion 
through a Rabi angle 9 determined by the energy and 
the length of the pulse. The excited state then decays 
via the three available channels. The 5D5/2 state will 
be maximally populated at a Rabi angle 9 equal to tt. 
By fitting values for the dark state probability at vari- 
ous values of the incident pulse energy with the expected 
form Pdark = B' sin^ (^/2), where 9 — a^/E, the branch- 
ing ratio B' for decays to 5D5/2 can be extracted. The 
fitting parameter corresponding to the excitation effi- 
ciency a, a function of beam waist and position of the 
ion in the intensity profile of the beam, was found to 
be 3.41 ± 0.07 nJ~^/^ and varied by no more than 10% 
from run to run. The fit parameter corresponding to the 



branching ratio B' proved robust to deviations in a at 
and beyond this level. After multiple runs at decreasing 
and increasing attenuation (corresponding to increasing 
or decreasing the energy per pulse), data were binned by 
pulse energy and fit to the two-parameter model above, 
the results of which are displayed in Fig. 4. Vertical error 
bars are statistical, corresponding to binomial standard 
error in the average of all runs. Systematic effects con- 
sidered in this measurement include the effect of adjacent 
pulses, which was shown through numerical integration 
of the optical Bloch equations to have little effect on the 
maximum of the dark state probability, and the uncer- 
tainty in shot-to-shot pulse energy, found to be approx- 
imately 10% which is indicated in the horizontal error 
bars. Near the maximum, this pulse energy variation 
contributes to a systematic shift in the measured value 
of B' of -0.7 %. The bright ion counts resulting from im- 
proper synchronization of the pulse picking AOM with 
the continuous pulse train of the Ti: Sapphire oscillator 
were counted during the course of data collection and 
found to occur an average of 4 ± 1% of the time and 
were removed from the data sets, leaving a 1% system- 
atic uncertainty. The same systematics as in the pre- 
vious measurement, namely the leakage of cooling light 
and collisions with background gas, again contribute ap- 
proximately 0.5% systematic uncertainty as previously 
mentioned. The combined effect of these systematic ef- 
fects leads to an overall systematic uncertainty of about 
2%. 
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FIG. 4: Data obtained from single pulse excitation of the 
6S1/2 — > 6P3/2 transition. Shelving probability represents 
the percentage of the runs where upon spontaneous decay, 
the ion was no longer fiuorescing, i. e. was in the 5D5/2 
state. The solid line is a two-parameter fit to the data us- 
ing Pdark = B' sin'^ {a.\/E/2^ , from which the value for the 
branching probability B' from 6P3/2 to 5D5/2 was extracted 
to be 0.215±0.0047. 

Combining the information from the weak shelving and 
the Rabi rotation experiments, one can calculate abso- 
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lute values for the three transition strengths. Data are 
summarized in Table (1) and pictorially represented with 
overall error bars in Fig. 5, and compared to previously 
measured values [T3J [THl [T71 (TB] as well as calculated ma- 
trix elements [TU [T21 [12] • The results of our measurement 
represent an approximate threefold improvement over 
previously achieved accuracy for two of the three transi- 
tions, with the notable exception of the 5D3/2 decay ratio 
whose error is approximately the same as the best pre- 
vious measurement. The transition probability Afi for 
each decay can be calculated using an experimentally de- 
termined lifetime tqp^^^ = 6.32(10)ns and the fact that 

Ti = {^f^fi^ ■ From these ratios, following the for- 
mal treatment of calculations of transition strengths [T4j 
115] ■ the radial matrix elements for the involved transi- 
tions can be calculated. We find (6P3/2I | -D^i ||6Si/2) — 
4.720±0.040, (6P3/2II II5D3/2) = 1.349 ±0.036 and 

(6P3/2II I>isi II5D5/2) = 3.945 ± 0.066 in units where 
eOo = 1. Our findings for the decay to the 6S1/2 state 
show agreement with the calculated values of [13] to 
within one standard deviation and are within the error 
bars of the [IH] and [T7] measurements. While our mea- 
surement of the 5D3/2 decay rate did not improve the un- 
certainty in the measured value for that particular decay, 
it is characteristically lower than all previously measured 
values, though not outside of quoted error bars. How- 
ever, all calculated values for that particular decay fall 
within the bounds of our measurement. Our value for the 
5D5/2 decay rate, which represents the most significant 
improvement in statistical uncertainty, sits at the lower 
bound of the [TB] and [17 measurements. 

Experimental confirmation of calculated values for 
these matrix elements is of particular use for several 
reasons. The parity-nonconserving electric dipole tran- 
sition amplitude for the transition 6S1/2 ^ 5D3/2 is 
given as a sum over matrix elements of both the dipole 
and weak charge Hamiltonian where the largest contribu- 
tions to the sums come from the 6P1/2 and 6P3/2 terms, 
approximately 90 and 8%, respectively [111 HO]. The 
(6P3/2II i)_Ei II5D3/2) and (6P3/2II -D_Ei II6S1/2) terms 
appear explicitly in this amplitude. Additional applica- 
tion of measured values for these matrix elements is of 
relevance to searches for drifts in the fine structure con- 
stant, both astronomically in anomalous quasar absorp- 
tion spectra |21j and in Earth-based optical frequency 
atomic clocks [55]. 

In conclusion, we have measured the branching ratios 
from the 6P3/2 state of ^^^Ba^ into all states allowed via 
dipole transitions using a single trapped ion excited with 
ultrafast optical pulses, and calculated transition prob- 
abilities and dipole matrix elements of the transitions. 
For two of the three transitions, we achieve a threefold 
improvement in overall uncertainty, while for the 5D3/2 
decay we are approximately at the same level of precision 
as the best previous measurements. The D state decays 
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FIG. 5: Comparison of presently measured transition rates for 
the three transitions from 6P3/2 to the final state noted at the 
left vertical axis to previously measured and calculated values. 
Error bars are one sigma confidence intervals on experimental 
data. Vertical ordering in lower two panels is the same as 
topmost panel. 



show a characteristically lower decay rate than previously 
measured and are closer to theoretical values. Ultrafast 
excitation of ions using single femtosecond pulses has par- 
ticular use for quantum computation when used with the 
137 isotope of Ba+. Coherent coupling between the spin 
states of the ion and the photon polarization state has 
been demonstrated in "^Cd^ [H] and ^'^^Yh'^ [T^ and 
its application to barium follows very closely the methods 
employed in the present work. 
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TABLE I: Calculation of transition probabilities and radial matrix elements from measured relative transition strengths from 
the 6P3/2 state in Ball and comparison to previous work. 



Transition 




Present work 




Previous Experiment 


Theory 




Rel. Strength 


D (a.u.) 


Af^ X 10'-' s"^ 


Afi X 10^ s-i 


D (a.u.) 










0.106±0.009[l3] 


4.6982[14] 


6P3/2 6S1/2 


0.756±0.046 


4.720±0.040 


0T196±0.0020 


0T17±0.004[18] 
0.118±0.00817. 


4.658[l5] 
4.6738[T^ 










0.00469±0.00029[T1] 


1.2836[Ti] 


6P3/2 5D3/2 


0.0290±0.0015 


1.349±0.036 


0.004589±0.00025 


0.0048±0.0005[IH] 
0.0048±0.0006[17] 


1.312[I1] 
1.2817[19] 










0.0377±0.0024[13] 


3.9876[14] 


6P3/2 5D5/2 


0.215±0.0064 


3.945±0.066 


0.03402±0.00115 


0.037±0.004ll8] 
0.037±0.004[T7] 


4.057[l5] 
3.8336[IS] 



